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Supersymmetry predicts the identity of Yukawa and gauge couplings in the QCD sector: qqg = qqg = qqg. We examine 
whether the qqg Yukawa coupling can be determined, by methods complementary to LHC, by analyzing squark-gluino 
final states at a TeV e+e — collider. 



1. Introduction 
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iy-j 'While in QCD non-Abelian gauge invariance leads to equal couplings for quarks and squarks to gluons, supersymmetry 
uarantees the identity of the Yukawa coupling between squarks, quarks and gluinos with the gauge coupling, cf. 
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9s {qqg) = g s (qqg) = g s {qqg) ■ (1) 

This identity is preserved if supersymmetry is broken by soft terms, i.e. gaugino/scalar masses, and bi/tri-linear 
couplings between scalar fields. The relation is crucial for the natural extension of the theory from the electroweak 
scale to the Planck-scale without introducing quadratic divergences, which are generated in the squark propagator, 
for instance, by the bosonic squark/gluon loops and which are canceled by the fermionic quark/gluino loop. 

Several methods can be exploited to test the identity ([T|) by measuring the magnitude of the Yukawa coupling qqg. 
At the LHC the production of squark pairs in quark-quark collisions which is mediated solely by the interchange of 
gluinos, provides a classical instrument for the measurement of the Yukawa coupling in SUSY-QCD Q|. In practice, 
however, an ensemble of auxiliary measurements of decay branching ratios is necessary, presumably requiring LC 
supplements, if the coupling should be determined in an [almost] model-independent way. 

A potential complement to this method is gluino emission in association with quark-squark final states in e + e~ 
collisions, cf. Refs. 0,0]: 

e + e" -> qqg (2) 
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1 INTRODUCTION 




Figure 1: Generic Born diagrams corresponding to the processes {Sp and jJP,RP which are proportional to the Yukawa coupling 
g s [box] and to the gauge coupling g s [dot], respectively. 



which requires a non-zero Yukawa coupling. This process is related by supersymmetry directly to the gauge processes 
of gluon radiation off squarks, 

e + e" -> qqg , (3) 

and standard gluon radiation off quarks, 

e + e" -> qqg . (4) 

Generic diagrams are displayed in Fig. [TJ [Note that the abbreviations qq.-, used for notational clarity, should 
anywhere be interpreted as the incoherent sum of particle-antiparticlc plus antiparticle-particle states of all 
flavors and L, R indices.] It can be anticipated without analyzing details that the process © will be useful 
in practice only if the squarks and gluinos are moderately light; the range, however, is quite compatible with 
electroweak precision analyses, cf. Ref. The process is suppressed by the [small] Yukawa coupling squared, 
the propagator transporting high (virtual) masses, and the 3-body phase-space involving two heavy particles; the 
suppression remains effective even significantly above the threshold. 1 Nonetheless, the production cross sections 
are in general expected still to be considerably larger than loop-mediated gluino-pair production in e + e _ collisions 0]. 

Since any method for measuring the SUSY-QCD Yukawa coupling is connected with ancillary problems, it is worth 
studying to what extent the process ([2]) could be exploited to approach a solution of this experimental challenge. 
Since the identity of Yukawa and gauge couplings is a central concept of supersymmetric theories, it is mandatory to 
study all opportunities which could shed light on this aspect. To reach solid conclusions, the theoretical analysis has 
been performed including SUSY-QCD corrections at next-to-leading order [NLO] which increase the complexity of 
the theoretical work enormously. In turn, this allows us to study the sensitivity of SUSY-QCD corrections in qq, qq 
final states to the Yukawa coupling. Preliminary results had been contributed to conference proceedings [lCJ . 



1 These difficulties contrast with measurements of the Yukawa couplings in the electroweak SU(2)xU(l) sector in which chargino, 
neutralino and selectron pair-production in polarized e+e — and e~e~ processes allow the determination of the electron-related Yukawa 
couplings at the per-cent down to the per-mill level §. 
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TTlq <C TTlg VTlq + nig < yJ~S < 2m g 

2m g < -y/s 


qq(2j + $) 
qq(2j + 9) 
qq(2j + $) 


qqg(3j + $) 
qqg{3j + $) 
qqg(3j + $) 


qqg(4j + $) 
qqg(4j + $) 


qqgg(6j + 9) 
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q<i(6j + m 


qqg(7j + If) 


qqg(6j + $) 
(qqg(6j + $)} 


qqgg{6j + $) 

(qqgg(6j + 9)) 
(qqgg(6j + #)) 



Table I: Kinematically accessible processes involving squarks and gluinos at different cm. energies for m,q < m g and rriq > m g . 
Processes are identified by the jet topology of the final states. The corresponding gauge and Yukawa couplings are denoted by 
a s — ffs/47r and a s — g^/Aiv, respectively. [Final states in (...) brackets in the second part of the table are generated primarily 
by squark decays to gluinos to which the power counting of the couplings does not apply. ] 



2. Squarks and gluinos in e + e collisions 

Before discussing the dependence of individual processes on the SUSY-QCD Yukawa coupling, an overview should 
be given on the size of cross sections which can naturally be expected for experimental analyses. In addition, 
the key points of the crucial techniques used in the higher-order calculations are summarized. For the reader's 
convenience, we include some earlier material from the literature as to provide a coherent presentation of the problems. 

The topology of the final states after the decays of the supersymmetric particles down to the LSP [assumed to be 
the lightest neutralino Xi m the present R-parity conserving set-up] depends strongly on whether squarks are lighter 
or heavier than gluinos. In the first scenario, squarks decay preferentially to charginos/neutralinos, while gluinos 
decay to squark-quark pairs, followed by the subsequent squark decays, 

m s < rrig : q -> q + X -> 3 + tf, etc 

. . „ ( 5 ) 

g -> q + q -> q + q + x -> 33 + tf, etc . 

In the second scenario squarks decay to gluinos which subsequently decay, mediated by virtual squarks, to 
charginos/neutralinos, 

niq> m g : q -> g + q -> q vilt + q + q -> q + q + q + x -» 3 33 + rf, etc 

. . „ ( 6 ) 
9 -> ?virt + q -> q + q + x ->■ 33 + 9, etc . 



In contrast to QCD jets, the jets in heavy squark/gluino decays are well separated in phase space, and clustering will 
occur only at random. A large fraction of the charginos/neutralinos x is m general light and, if not the LSP, they de- 
cay to leptons [with taus generating only slim jets with up to three hadrons which we keep separate in the jet counting.] 



A survey of all squark and gluino production processes relevant for investigating the strong gauge and Yukawa 
couplings is given in Tab. [T] which indicates the threshold energies and the jet topologies for the individual continuum 
processes in which q or g are not mutual decay products [QCD coupling a s = gl/Air; Yukawa coupling a s = g^/Air]. 
At the time when the measurement of the Yukawa coupling can be performed, the masses of the supersymmetric 
particles will be determined very precisely, c/. Ref. and it will be evident which part of the table is relevant. 



It may be noted that final states of events generated by the decay of on-shell squarks to quarks and the lightest 
neutralino as LSP, q — > q + Xii can be reconstructed up to a two- fold ambiguity, an important identification 
tool when the squark and neutralino masses are known, c.f. Ref. As the qq cm. frame [denoted by *] 
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coincides with the laboratory frame, or, can be reconstructed after gluon emission, the qq axis in that frame is 
given by one of the two intersections of the cones centered at the quark axes with opening angles determined 
from m| = m~ — 2E*E~(1 — /3|cos0|). The false solutions for the axes give rise to essentially flat background 
angular-correlations. 




[GeV] [GeV] 

Figure 2: Relevant tree level cross sections for the production of squarks and gluinos. In Fig. (a) the gluino is assumed to be 
heavier than the squarks, according to the SPSla scenario, and in Fig.(b) the reversed mass scenario is considered, where the 
gluino is lighter (trig = 530 GeV,) than the squarks (degenerate squark masses uiq — 600 GeV,). 

Leaving out rare processes of second order in the Yukawa coupling, the cross sections for squarks lighter than 
the gluino are displayed in Fig. [DJa), summed over five flavors. For clear identification a minimal energy cut 
Ex > 100 GeV is applied to the gluon jet (details are given in the following subsections). In Fig. [5]Ja) produc- 
tion cross sections are shown for squark and gluino masses adopted from the standard SPSla scenario [131 ], i.e. 
nig = 600 GeV and a common squark mass of rriq = 530 GeV, [which is compatible with the mass limits allowed 
by the searches for supersymmetric particles at the Tevatron [14|]. The dominant cross section of the qq final state 
rises steeply above the threshold proportional to the third power of the velocity of the squarks. Since the gluino is 
heavier than all squarks in the SPSla scenario, the production threshold for the qqg final state is higher than the 
qq threshold. The jet topologies in the final states are characteristic for the individual processes and jet counting 
can be used as a powerful discriminant. [Of course, this Born approach must be refined by properly controlling QCD 
showers in experimental analyses.] 

The relevant cross sections for the scenario in which the squarks are heavier than the gluinos, are shown 
in Fig. [Hb). In a merely ad hoc procedure, as any details beyond SUSY-QCD are irrelevant for the present 
discussion, the mass values of squarks and gluino in SPSla are just reversed. [Such a scenario is actually close 
to SPS4.] Note that above the cm. energy of twice the squark mass the process qqg proceeds primarily through 
squark pair production with subsequent decay of a squark into a gluino. The energy window for qqg production sui 
generis is narrow and the cross section of the process, due to phase space suppression near the threshold, is rather small. 

However, an additional opportunity for measuring the Yukawa coupling in the second scenario is provided by the 
squark decays to gluinos, Ref. [3]. The threshold behavior of the squark-excitation cross section in e + e _ collisions 
depends sensitively on the total decay width r(g) tot . This method of measuring the total width was discussed in 
Ref. 0] for slepton production but it can be transferred to squarks in the same way; in fact, the impact of the total 
width on the production cross section for squarks is significantly bigger as a result of the increased size of the width. 
On the other hand, the electroweak continuum cross section for squark-pair production with squarks decaying to 
gluinos is proportional to the branching ratio BR(q — > qg) squared. By combining the two threshold and continuum 
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measurements of total width and branching ratio, respectively, 



r(<? -> qg) = BR{q -> g<?) x r(g) to t = 1 1 - m 2 

the SUSY-QCD Yukawa coupling can be determined from the partial width. 



(7) 



Since in the analysis of the individual channels higher-order SUSY-QCD corrections will be included, the technical 
set-up should be summarized globally. To take advantage of our extensive infra-structure, the dimensional regulariza- 
tion has been carried out in the DREG scheme. To comply with the SUSY-QCD Ward identities, which are violated 
by the mismatch of fcrmionic and bosonic degrees of freedom in this scheme, counter terms must be added which 
however can be mapped onto a modified relation between gauge and Yukawa couplings 



15 



3: 



9s = 9b 



1 



a a (2 



9s 



(8) 



The masses of the heavy particles, squarks and gluinos, are introduced as pole masses. For the renormalization of 
the couplings the MS scheme is adopted. The "experimental" QCD gauge coupling a s 5 \Q 2 ) is introduced for five 
massless quarks and all heavy particles are decoupled. From this coupling, conventionally fixed at the energy Q = Mz, 
Rcf. H3, the appropriate QCD gauge coupling for the SUSY system, a s (Q 2 ), above the heavy particle thresholds can 
be derived as 



a s (Q 2 



a 



(5) 



(Q 2 




(9) 



with the number of flavors Nf = 6 and colors Nc = 3. This procedure guarantees that the coupling a s (Q 2 ) smoothly 
approaches the asymptotic form when Q 2 ^> M 2 for all masses M of the supersymmetric theory. In the calculations 
presented here we identify the renormalization scale with the cm. energy. 



After this introductory discussion of the individual channels and their dependence on the SUSY-QCD Yukawa 
coupling, detailed evaluations are presented in the following subsections. 



2.1. Final-state supersymmetric particles 

2.1 .1 . Squark pairs qq 

It has been shown in Figs. [2] that the cross section for squark-pair production 

e+e- -> qq* (10) 

is maximal in the group of channels suited for the measurement of the SUSY-QCD Yukawa coupling. The cross 
section, in Born approximation, 

7TCX 2 

VLo(e + e~ -> q L/R q* L/R ) = —C L/R f3 3 , (11) 

rises as the third power of the q velocity /3 = (1 — 4m'~/s) 1 / 2 above the threshold as demanded for P-wave production 
0. The L/R couplings, 

Cl/r = Q 2 q + (c e v 2 + c\ 2 ){c q v ± c q A ) 2 X % - 2Q q c L v (c q v ± c\)xz (12) 

are given by the electric and isospin charges of the squarks, Cy = (// — 2QfS 2 v )/s2\v and = I^/s2W, where 
S\y = sm 2 0\y is the electroweak mixing parameter, s 2 w — sin 28yy, and xz = s /( s ~ m %) the scaled Z-propagator. 
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vertex correction, mid: gluon bremsstrahlung; right: gluino vertex correction. 



The production cross section for squark pairs in next-to-leading order (NLO) includes the gluon vertex correction 
and gluon bremsstrahlung, as well as the gluino vertex correction, as depicted in Fig. [31 [Final states qqg are assumed 
to be separated explicitly.] The size of these corrections, 



qq 



qq*g) = cr(e + e _ 



qq*) 



I Born 



° Ls A vert +real , _ ^ s a vert 
3 7T S lu ° n 3 7T SluIn 



(13) 



depends on the masses of the squarks and gluinos, cf. Ref. 18]. The two parts contributing to the corrections 
are displayed, for the set of SPSla parameters defined earlier, in Fig. 3J The gluon vertex contribution is formally 
divergent at the threshold, i.e. linearly in the inverse velocity of the squarks, but this Coulomb singularity is 
regularized by the non-zero squark width 19j. [In toto, the Born threshold suppression ~ f3 3 is reduced to ~ f3 2 .] 
Compared to the Born cross section, O(50 fb), the corrections, in particular the gluino corrections involving the 
Yukawa coupling, are of order femtobarn. With an integrated luminosity ~ 1 ab _1 , these corrections can clearly 
manifest themselves beyond statistical fluctuations. 
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Figure 4: The NLO SUSY-QCD (a s and & s ) corrections for sguark pair production. Separately shown are the sum of virtual 
and real gluon corrections, (4:a B /3iv)A v J^ t * real , as well as the gluino vertex correction, (4a s /3ir)Agf^ no . Parameters correspond 
to the reference point SPSla. 



2.1 .2. Squark-pair plus gluon final state qqg 

A subset of final states in I2.1.1I is generated by the process 



e + e 



— * 



(14) 



in which the gluon jet is isolated explicitly. This process measures the QCD coupling of gluons to squarks. The 
differential cross section in the scaled Dalitz variables Xq = 2Eq/ ' \fs, etc, was calculated already some time ago (2p| : 



dcr(e + e -> q L / R q"* L / R ff) 



4.s 



-C 



L/R 



16 



16 



2 1 



16 Mi*y|* 



3 (1 



8(2/i|-l)(y| + y|, -y 2 g ) 
3 



(l-Xg){l-X r ) 



(15) 
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To simplify the notation, the abbreviations y q — 2\p q \/y/s, /Ug = rriq/y/s, etc, have been introduced. If the squarks 
are lighter than the gluinos, as realized in SPSla, the final state consists of three well-separated jets and a pair of 
ncutralinos /charginos. 

In the numerical example, Fig. [6j we assume a minimal energy to be carried by the jet(s) emitted in addition to 
the squark pair. The minimal energy Ex > 100 GeV required for X in the final state qqX cuts out the infrared 
divergence in the gluon energy for X = g in Born approximation. The dashed line in Fig. [6] corresponds to the Born 
approximation. The solid line adds the complete set of SUSY-QCD radiative corrections. These corrections include 
vertex and box diagrams as well as additional gluon-pair and light quark-antiquark final states. For these final states, 
X = gg and qq, the same cut-off in the energy (Ex > 100 GeV) is applied as before. Example diagrams for the 
radiative corrections are shown in Fig. El The semi-inclusive definition integrates out both the infrared and collinear 
singularities in the NLO corrected cross section. 

In Born approximation the cross section depends strongly on the renormalization scale introduced solely by the 



q A A 




Figure 5: Example diagrams for SUSY-QCD corrections to squark pair plus gluon production, left: gluino vertex correction, 
mid: box diagram; right: double gluon bremsstrahlung. 



gauge coupling. The dependence is significantly reduced if the higher order corrections are included. For example, 
varying the renormalization scale /i between hsfs < /J, < 5y/s around a cm. energy of t/s = 3 TeV, the cross section 
varies by Aa/a ~ 14% in Born approximation, while the variation is damped to 4% in NLO approximation. 

In the peak region we find a sizeable cross section of several femtobarn. The radiative corrections enhance the cross 
section at a level of 20%. 
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7i [GeV] 

Figure 6: Total cross section of the process e + e~ — > qqg for the cm. energy y/s at leading order [dashed line] and next-to-leading 
order [full line]. The mass parameters are adopted from SPSla. A minimal energy Ex > 100 GeV is assumed for single jets g 
and jet pairs gg, qq emitted off the squarks. 
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2.2. The golden channel: Squark-gluino final states qqg 



The golden channel for measuring the SUSY-QCD qqg Yukawa coupling g s directly is the process 



e + e — > qq*g + c.c. 



(16) 



As evident from Fig. [H the size of the cross section is governed by the quadratic dependence on the Yukawa coupling 
g 2 = Aira s : 



d 2 a 2a 2 a< 



3.s 



C 



L/R 



2/1 2 i 2 \ 



(1 )2 [3(1 - as,)(l ~ asg) + 3/xf - 3/xf. - a*(/i§ - - ^(1 - s,)(t£ - 2/1 + $)] (17) 



1 



5 "t"q' ^q\^g ^q* 

.2 o,,2 \(„3 „,2 , „,2 



(l-x g )(l-a; $ .) 
in the standard notation introduced before. 



qyr - 9 ( x a- + 2 4 - M-m -y§ + v?) 



The Born cross sections for both the cases vriq < m~ g and rriq > m,g are presented in Figs. [21 The processes give rise 
to 4, respectively 6 jets in the final state [at Born level], with the 4-jet final state characteristic for qqg production in 
the mass range rriq < rrig. 

Examples of NLO diagrams contributing to the SUSY-QCD corrections to qqg final states are depicted in Fig. [7] 




Figure 7: Example diagrams for SUSY-QCD corrections to squark gluino production, left: gluon vertex correction, mid: box 
diagram; right: real gluino emission. 

Vertex and box diagrams are built up by gluinos, gluons, squarks and quarks. To order a s , 4-parton final states 
are generated only by the emission of an additional gluon. Since the radiatively corrected cross section is defined 
inclusively infrared and collinear singularities are integrated out. 

The corrected cross sections are displayed in Figs. [5] for squarks lighter and heavier than gluinos [SPSla and 
SPSla reversed], respectively. Again, the strong dependence on the renormalization scale in Born approximation 
is damped once the higher order corrections are included properly: If the renormalization scale \i varies between 
\\fs < [i < 5y/s, for ^/s = 3 TeV the cross section for rriq < nig varies by Aa/a ~ 14% in Born approximation, while 
the variation is damped to 4% in NLO approximation. 



The cross sections for rriq < ni g are enhanced by the radiative corrections to a size of order 2 fb which can be 
considered sufficient for experimental analysis if integrated luminosities ~ 1 ab _1 can be achieved. In the mass range 
rriq > nig the size of the cross section is rather small, however. 
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2.2 Squark-gluino final states qqg 
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Figure 8: (a) Total cross section for squark-gluino production at LO and NLO for SPSla mass parameters, i.e. trig < m g ; 
(b) for reversed SPSla mass parameters, i.e. mq > nig. [Above a cm. energy of twice the squark mass, squark-pair production 
followed by squark decays to quarks and gluinos, is by far dominant.] 



2.3. Final states without supersymmetric particles 

2.3.1 . Quark pairs qq 

The basic quark-pair production process, defined in the absence of any supersymmetric particle in the final state, 

e+e- -► qq (18) 

scales asymptotically with the energy squared: 

Ana 2 



a(e + e -> qq)\ 



f q + (c e 2 



)Xz - 2Q q c e v c q v xz 



(19) 




Figure 9: Examples for diagrams of QCD and SUSY-QCD corrections to the production of quark pairs; left: gluon vertex 
correction, mid: gluon bremsstrahlung; right: squark/ gluino vertex correction. 

The Born cross section is modified, to NLO, by gluon exchanges in the vertex and by gluon bremsstrahlung, as well 
as squark/gluino corrections to the vertex, Fig. [51 



CTNLo(e + e -> qq) = a(e + e -> qq) 



Born 



1 



_ a vcrt+rcal 
3 7T glu ° n 



4 <Xs_ a vert 

3 n ^squ/gluino 



(20) 



The gluon correction to NLO is well known, A^*^ rcal = 3/4, while the genuine triangular squark/gluino correction 



is more involved 2l|. For the set of SPSla and reversed SPSla mass parameters introduced before, the pure QCD 
corrections are in general leading, except for very high energies. The individual gluon and SUSY corrections are 
presented in Fig.[TUJ [The two contributions are separately gauge invariant.] The Born cross section drops from 2.2 pb 
down to 61 fb in the range of cm. energies between 0.5 and 3 TeV. The genuine SUSY-QCD corrections involving 
the Yukawa coupling of squarks to quarks and gluinos are moderately large and of half the size of the standard QCD 
corrections. In the peak region the radiative SUSY corrections are of the order of 1.5%; the corrections flip sign at 
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Figure 10: NLO QCD, (4q s /3tt) A v g ^ + real , as well as NLO SUSY-QCD, (4a,/3ir)A^ 9lttjno) corrections o/ the total cross 
section for quark pair production in the reference points SPSla and reversed SPSla. 



an energy of ~ 3 TeV for the two reference sets of SUSY masses. For rriq = nig = M the correction can be cast in 
the compact form 



A vert JJJW 
squ/gluino 4 4 I 1 -4- Q 



In'' 



1-/9 
1 + /3 



where /3 = (1 — 4M 2 /s) 1 / 2 . The SUSY correction approaches 

= -1/4 ln(s/M 2 ) 



a vert 
squ/gluino 



(21) 



(22) 



for asymptotic energies. The logarithmic behavior in the energy is consistent with the Lee-Nauenberg theorem since 
gluino radiation is not included in the final state; in fact, adding gluino radiation cancels the logarithm in the 
asymptotic region. 

2.3.2. Quark Pairs and Gluons qqg 

If additional gluon jets are isolated in the final state, 



qqg, 



the analysis becomes much more complex. Starting from the Dalitz plot density for gluon radiation 22j, 



d 2 a{e+e- 



qqg) 



a{e + e -»■ qq)\ 



2a, 



Born 



(23) 



(24) 



Born 3lT (i-^Xl-Zf) ' 

with the scaled energies x q = 2E q /^/s etc, the jet events are defined for a cut in the invariant masses of the parton 
pairs yij — (pi + pj) 2 /s > 0.05. The radiative corrections may be decomposed in the form 



0NLo(e + e" 



qqg) 



a(e + e 



qqg) 



1 



. vert+rcal 



; a vert 

squ/gluino 



(25) 



As generally anticipated for corrections to jet cross-sections, the gluon corrections, c.f. Ref. 23 1 



—A 



vert+real 
gluon 



(26) 



a s B 

7T s '""" 7T 2 A 

are by far dominant with B/2A ~ 10 for y tj = 0.05, i.e. ^A^*+ real ~ 0.3. The genuine SUSY-QCD corrections are 
defined as any corrections involving propagators of supersymmetric particles, some examples are shown in Fig. HU 
Such diagrams come with two powers of the quark-squark-gluino Yukawa coupling at next-to-leading order. The 
genuine SUSY contribution ^ ^squ/ gluino ^ s displayed in Fig. 1121 The two mass hierarchies rriq < rrig and rriq > rrig 
are adopted again from the reference point SPSla and the reversed SPSla point. In the respective peak regions SUSY 
corrections of about 1.5% are predicted, significantly smaller than the standard QCD corrections, as anticipated. 
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Figure 11: Examples for diagrams of SUSY-QCD corrections to the production of quark pairs accompanied by a gluon; left: 
squark/gluino vertex correction, mid and right: squark/gluino box correction. 
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Figure 12: NLO SUSY-QCD corrections, (a s / ^)^ S q U/ /gi U i no , to the total cross section for qqg production for the mass hierarchies 
of SPSla, solid line, and SPSla reversed, dashed line. 



3. Assessment 



In the preceding section the theoretical base has been studied for future measurements of the Yukawa coupling 
between quarks, squarks and gluinos in high-energy e + e~ collisions. Supersymmetry predicts this coupling to be 
identical with the standard QCD gauge coupling between quarks and gluons, and squarks and gluons. Various 
channels have been investigated which can be explored to test this fundamental identity between the couplings, 
either indirectly through virtual corrections to quark and squark pair production in e + e~ collisions, or directly in 
comparing the cross section for the golden channel qqg with the standard gluon radiation processes qqg and qqg. 
While the golden channel measures the qqg Yukawa coupling g Sl the radiation processes qqg and qqg determine the 
QCD gauge coupling g s in the squark sector and the standard quark sector for comparison. 

The central results are summarized in Tab. HT1 and Figs.[l3l In the table the SUSY parameters are adopted from 
the reference point SPSla with squark and gluino masses of 530 GeV and 600 GeV, respectively. Since by the time 
of measurements of the couplings the squark and gluino masses will be determined accurately, special values of the 
e + e~ cm. energy can be chosen such as to maximize the sensitivity for the measurements of the Yukawa and gauge 
couplings. A standard value j£ = 1 ab _1 is taken for the integrated luminosity. The errors quoted in the table are 
purely statistical, reflecting the minimal errors that can be expected in any future experimental analysis. 



The left-most figure in the panel Figs. [13] displays the contours of the statistical 1-sigma errors Ad s /d s for the 
standard qq channel in the [m g -,r«j] mass plane for ILC, irrespective of any specific model except SUSY-QCD. For 
each g, g mass pair the collider energy is chosen such as to provide maximal sensitivity to the measurement of a s ■ 
The Yukawa coupling d s can be determined in this channel only indirectly through virtual SUSY corrections. Such 
analyses may provide a first glimpse of the Yukawa coupling. 
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Table II: 1-sigma statistical errors of the Yukawa a a and gauge ct s couplings determined from quark- s quark- gluino final states 
and SUSY radiative corrections in various channels based on SPSla SUSY parameters. The upper part of the table with the 
final states qq, qq and qqg describes the measurement of & s , and the lower part, with qqg and qqg, the measurement of a 3 . 
The cm. energy y/s{max\ is chosen, below 3 TeV, such that the sensitivity to the measurement of the couplings is maximized, 
except for the gauge coupling in qqg where the scale is adjusted to qqg. 



e+e~ ->• qq, y/s < 1 TeV e+e~ -> qqg, y/s < 1 TeV e+e~ -» qqg, y/s < 3 TeV 
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Figure 13: Contours of the statistical 1-sigma errors of the Yukawa coupling Aa s /a s in the [mq,mg] mass plane, Left: for 
the indirect channel e + e~ — ► qq at ILC energies; Mid, Right: for the direct golden channel e + e~~ — > qqg, corresponding to a 
maximal ILC cm. energy of 1 TeV, and to 3 TeV for CLIC, respectively. The blind wedge can be covered by analyzing the 
partial widths. 



The figures in the middle and on the right of the panel summarize statistical analyses of the golden qqg channel 
on the same basis. In this channel the Yukawa coupling can be measured directly since the Born cross section is 
proportional to a s . The results are presented for ILC cm. energies up to a maximum of 1 TeV, and for CLIC up 
to 3 TeV. As before, the energies are chosen at the point of maximal sensitivity for the Yukawa coupling for any 
mass pair [m,q,mg]. For m g < rriq the collider energy is restricted to the range below the qq threshold as to avoid 
contamination of the qqg channel due to on-shell q — > qg decays. This restriction generates a blind wedge near the 
threshold which however can be covered in analyses of the partial widths as discussed before. 

The estimate of the statistical error Ag s /g s ~ 1% for the measurement of the Yukawa coupling in this report is 
only a first step. The step is promising, nevertheless, as the expected error is sufficiently small. Even if an order of 
magnitude of events is lost due to experimental cuts and efficiencies in controlling, for example, the chargino and 
neutralino modes of squark decays [and cascading gluinos], the resulting final error is still at a level of 5 to 10%. This 
size compares reasonably well with expectations for LHC experimental analyses in Ref. [4| in which squark/gluino 
decays are included explicitly. Thus, LHC and a TeV lepton collider are complementary instruments for the involved 
measurements of the SUSY-QCD Yukawa and gauge couplings. In particular, since the relevant cross sections drop 
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with rising squark and gluino masses quadratically and faster 2 , mutual complementarity of the two methods raises 
prospects of testing the identity of Yukawa and gauge couplings in super-QCD considerably. 
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